Aging in the glass phase of a 2D random periodic elastic system 
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Using RG we investigate the non-equilibrium relaxation of the (Cardy-Ostlund) 2D random Sine- 
Gordon model, which describes pinned arrays of lines. Its statics exhibits a marginal (6 = 0) glass 
phase for T < T g described by a line of fixed points. We obtain the universal scaling functions 
for two-time dynamical response and correlations near T g for various initial conditions, as well 
as the autocorrelation exponent. The fluctuation dissipation ratio is found to be non-trivial and 
continuously dependent on T. 

PACS numbers: 



Tremendous progress has been achieved in recent years 
in the detailed understanding of the off equilibrium re- 
laxational dynamics (e.g. coarsening) of pure systems. 
For instance the scaling forms of response and corre- 
lations are well characterized and several autocorrela- 
tion and persistence exponents have been computed us- 
ing powerful methods By contrast, analytical studies 
of aging properties of random and complex systems have 
been mostly achieved within mean field theory Q • These 
works have unveiled the existence of non trivial fluctu- 
ation dissipation ratios (FDR), X. These generalize to 
non-equilibrium situations the fluctuation dissipation re- 
lation between (disorder averaged) integrated response 
and correlation 2] and can be interpreted as an effective 
temperature T c g — T/X These were later inves- 
tigated beyond mean field for pure systems. Although 
they appear to be trivial for some pure domain growth 
processes Q , they were later found to be novel universal 
quantities at pure critical points 0, Q . 

The description of glassy states beyond mean field 
needs to take into account very broad distribution of re- 
laxation times originating from rare sample to sample 
fluctuations and distribution of barrier heights Q. In 
this context it is not obvious how (or whether) the off 
equilibrium dynamical scaling valid for pure systems ex- 
tends to this case. Also it is not obvious that various def- 
initions of FDR (using different observables, or different 
configurational averages) are equivalent and meaningful 
Hi H E! ■ Some results were obtained for the random 
mass Ising critical point, via RG. The dynamical (z) and 
autocorrelation (6c) exponent were computed A 
definition of the FDR (for the zero mode) was found to 
be non trivial 12] and obtained to one loop. 

Whether such results hold within a glass phase, away 
from a phase transition, is still an outstanding question. 
A prototype system which exhibits a glass phase, while 
analytically more tractable, is the disordered elastic man- 
ifold in a random potential 0. It describes a variety of 
physical situations, both in or out of equilibrium, rang- 
ing from visualization of pinned domain wall relaxation 
in magnets [lil non equilibrium transport in elec- 



tronic glasses 0] to mesoscopic fluctuations in vortex 
physics 0]. Here we investigate the case of a periodic 
manifold, a case which can be mapped 0, 0] onto the 
so-called Cardy-Ostlund (CO), 2D- random Sine Gordon 
model, defined by the hamiltonian: 

H[cb] = J d 2 x(\(V x <t>(x)f - Re(e(x)e^W) (1) 

where 4>(x) e] — 00, +00 [ is an XY phase (excluded vor- 
tices), £(x) is a quenched gaussian random field, i.e. a 
complex variable, with random phase and £(x)£(x')* = 
g5^ 2 \x— x'). The statics of Q has been extensively stud- 
ied analytically and numerically, in the context of planar 
flux line arrays (with displacement field it = a<f>/(2'K) 
and mean spacing a) and solid on solid models with dis- 
ordered substrates. It is known to exhibit a glass phase 
0, below T g (= 4tt ) described by a line of fixed 
points perturbatively controlled by the small parameter 
t = (T g — T)/T g (see however 23]). By contrast, its 
(non driven) non equilibrium dynamics has only been 
studied within gaussian variational approximation GVA 
[2^|. known already in statics to yield unreliable results 
for some observables 0|. There is thus the need for a 
controlled study via RG. 

In this paper we study the relaxation dynamics for 
the CO model in the glass phase T < T g starting from 
various initial conditions at tj = (e.g. with the same 
correlations as the pure model at equilibrium at temper- 
ature T"). We compute, to lowest order in r using RG 
along the fixed line, the two time (f < t) response 1Z, as 
well as the connected C and, respectively, off-diagonal C 
correlations. These are found to be characterized by (In- 
dependent) universal scaling functions of q z (t — t') and 
t/t', q being the wavevector. We find that an autocor- 
relation exponent can be consistently defined, i.e. with 
@c = ®Ri om y from the connected correlation (while the 
decay of C is too slow). Similarly, the associated FDR, 
found to be a scaling function, reaches at large time sep- 
aration (t/t' — > 00) a non trivial (and to one loop, q- 
independent) limit X^ which is non zero only for the 
connected correlation. This sheds light on relevant defi- 
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nition of the FDR in a glass phase. 

Let us start by a short discussion of the statics. The 
2-point correlation function exhibits anomalous growth 



{((p(x) - 0(O)) 2 ) = A{t) \n\x) + OQnx) 



(2) 



The RG to one loop [26j yields the universal 1271 A(t) = 
2r 2 + 0(t 3 ). Numerical simulations near T c [H and at 
T = agree qualitatively, and yield A(t — 1) = 02 = 
0.57 H3 and A(t = 1) = 2{2n) 2 B = 0.51 Hfl. A re- 
cent work j3]| claims an exact solution for the 2-point 
correlation using a conjectured mapping onto disordered 
free fermion model. The result, translated in the present 
variables would yield A(t) = 2r 2 (l — r 2 ). This however 
is clearly incompatible with numerics, showing that there 
is more (non perturbative ?) physics to be understood 
even in the statics of this model. A more direct signa- 
ture of the glassy nature of the phase is the sample to 
sample susceptibility fluctuations j2^ which are found 
(within one loop RG) to be universal and ~ 0{t) along 
the fixed line. The free energy exponent in this glass 
phase is 9 = 0, indicating a free energy landscape with 
logarithmic roughness. This is consistent with the find- 
ings in the equilibrium dynamics: an anomalous diffusion 
exponent [32|, continuously varying along the fixed line 
has been computed in RG which indicates a logarithmic 
growth of energy barriers with scale. These properties 
are caracteristic of a marginal glass (by contrast with 
the case 9 > described by a T = fixed point), ft is 
reminiscent of a related and simpler case of a vortex in a 
random gauge XY model where a freezing transition (at 
z c = 4) was found along a line of fixed points 33] . 

The relaxational dynamics of the CO model JU is de- 
scribed by a Langevin equation 



SH[4(x,t)} 



COM) 



(3) 



6<j)(x, t) 

where the thermal noise ((x,t) is such that (((x,t)) = 0, 
(C(x,t)C(x' ,t')) = 2r]T8( 2 \x-x')8{t-t') where T < T g 
is the temperature and rj the friction (77 = 1 in the follow- 
ing). The system at initial time t{ (= 0) is prepared in an 
equilibrium state of JJ| without disorder at temperature 
T = eT, [<j) q ,t=o<l>-q,t=o]i = T'q~ 2 , 4> q j being the Fourier 
transform, w.r.t. space coordinates, of the field <j){x,t). 
Since the disorder is irrelevant above T g , this choice of 
initial condition for T' > T g describes a quench from a 
high temperature phase to the glass phase (T < T„). A 
quench to high temperature (T > T g ) is studied in [34j. 

We will focus on the correlation C q t , and the connected 
(w.r.t. the thermal fluctuations) correlation C q t , 



C q - 
C q - 



[(<t>q t <t>-qt>) 



qtp-qt> )\i 



[{4>qt){4>-qt')]i 



and the response lZ g t , to a small external field /_ 



qt> 



K q 



~ S(<t>qtY 



t > t' 



(4) 



(5) 



where — , (..) and [..].; denote respectively averages w.r.t. 
disorder, thermal fluctuations and initial conditions. We 
focus on the FDR X q t , associated to the observable <j) 0] : 



(6) 



defined such that X q t , — 1 at equilibrium, i.e. when 
response and correlations depend only ont — t' [35j ] . 

The dynamics (j2J) of the CO mod el ffl is then studied 
using the standard MSR formalism |36tl37l] , using the Ito 
prescription. The correlations and response JSJ are 
then obtained from a dynamical (disorder averaged) gen- 
erating functional or, equivalently, as functional deriva- 
tives of the corresponding dynamical effective action T. 
This functional can be perturbatively computed [3^ us- 
ing the Exact RG equation associated to the multi-local 
operators expansion introduced in |39j . It allows to han- 
dle arbitrary cutoff functions c(g 2 /2Ap) and check univer- 
sality, independence w.r.t. c{x) and the ultraviolet scale 
Ao. It describes the evolution of Y when an additional 
infrared cut-off A; is lowered from Ao to its final value 
A; — > where a fixed point of order 0(t) is reached. In 
this limit, one obtains lZ q t , and C q t , (for t > t') from 



= 5{t-t') 
(7) 



dii / dt 2 K q ttl D c tit2 K q t , t 



Jti Jti Jti 

(8) 

where the self energy T, tl t 2 and the noise-disorder kernel 
D% t are directly obtained from Y at the fixed point [ioj ]. 
One finds (for details see |38ll4l|): 



E«/ - -2 £ - 1 e^ri?-7 e-5«/+</) (9) 
D\ t , = TqT-'e^re-i^+^il - e ~ c ^) (10) 

where je is the Euler constant, B^) = c[f* + c[f}, — 
2C^), C[f) = 7(i + t') - j(\t - t'\) is the bare Dirichlet 
propagator at coinciding points, B^f = e(2j(t + t') — 
7(i) — 7(t')), arising from the average over the initial 
condition. We have defined j(t) = J a c(a) ln(Ajjt + §), 
using the parameterization c(x) = J c{a)e~ ax for the 
cutoff function, and denote J2f t 7° = 9{t—t') f c(a)(A%(t— 
t') + |) _1 the bare response at coinciding points. Up to 
a boundary term, the correlation C q t , satisfies (jSJ setting 

-C id) 

the e «' term to zero in the above expression for D^ t , . 

Studying Eq. JJJ (setting ti — 0) , in the scaling regime 
q/A < 1, A 2 ,*, A^t' > 1 and keeping: 



v~q z (t-t') , u 



t 

T 



(11) 
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fixed, one finds a solution consistent with the scal- 
ing form: 



- F R (q*(t-t'),t/t') (12) 



similar to the form obtained for the response for critical 
systems ^HU^. Here the two exponents z and Or are 
identified from the logarithmic singularities of the scal- 
ing function, respectively at u — > 1 and u — > oo. We find 
that z identifies with the equilibrium dynamical exponent 
z — 2 = 2e 7E r + 0(r 2 ) and is thus, as expected indepen- 
dent of initial conditions under study. The exponent Or 
in (|12|l . characteristic of long time non-equilibrium relax- 
ation is also independent of e: 



r 2 ) 



(13) 



One also finds, with this choice of Or, that 
liniu—Kx, Fr(v, u) — Froo^). The full scaling function 
Fr(v, u) however depends on the initial conditions and is 
universal (up to a single overall non-universal length scale 
q — > Xq). Its explicit expression, given in |4lj . contains 
both non-equilibrium and equilibrium (u — > 1) regimes. 
We have checked to this order that it naturally splits 
into Fr(v,u) = F R q (v) + F^ oneq (u,u). Here we give the 
non-equilibrium scaling function only in the large time 
separation limit u — > oo: 



dt 2 / ^i(e t2 - 41 - l)£ tlt3 +<7 
(i Jo 

1— e \ 



1 



J tlt2 



t 2 +*1 



(i 2 -ti) 2 \\2^hh 



-1 (14) 



takes a scaling form C% 



with F^,{v,u) = e 



Tq- 2 F?(q 2 (t - t'),t/t') 

Sr, z = 2. The FDR 

X^ t , = (1 + e^ 2 "^)^ 1 interpolates between 1 and 1/2 
as u — t/t' increases from 1 to oo (for q ^ 0), while 
= \ the "random walk" value. In presence of disor- 
der one can solve (JSJ perturbatively in r using the above 
solution for the response. One obtains that C£,, in the 
scaling regime, is consistent with the scaling form 



C q 



The calculation [41| shows that q — Or, yielding an 



autocorrelation exponent \; 



= -40c 

0(r 2 ) 



1 



(17) 



, uF c {v, u) 



This value of 0^ is also such that lim„ 
F Coo {v). We have obtained the complete expression of 
F c (v,u) |4l| from which we can extract the large u be- 
haviour. We then discover the relation, valid for any e 

F Coo( v ) = (2 + ^ e t)vF Roo {v) + 0(t 2 ) (18) 

Having determined both response and correlation we 
now obtain, in the scaling regime, the FDR ©. It is also 
characterized by the universal scaling function 



= F x (q*(t-t'),t/t') 



(19) 



which also depends on e and has a complicated form [41| . 
In the limit of large time separation i.e. q z (t — t') fixed, 
t/t 1 — > oo, using l|14U18|) it simplifies into: 



where the constant a — J^ 3 dt% J Q dtil^ntz is a monotonic 
decreasing function of e. The large v behaviour is a power 
law. The form (|14f> is convergent due to the substraction 
of the pole at t 2 =t\. The logarithmic divergence asso- 
ciated to this pole yields a non trivial z exponent. The 
substracted piece precisely gives the equilibrium scaling 
function for the response F R q (v), up to 0(t 2 ): 

F R q (v) = e- v + re^ {{v - 1) Ei (v) e - v + e~ v - 1) (15) 

where Ei (v) is the exponential integral. The same result 
is also obtained if tj, is taken to be ij = — oo from the 
outset (e.g. at large but fixed system size), showing that 
the nonequilibrium regime of the scaling function merges 
smoothly with the equilibrium one |43j . 

We now turn to correlation functions. To obtain 
the equilibrium correlation one can simply use FDT 
which holds in this regime (i.e. X^ t , = 1) and Cj? t °/ q = 
Tq- 2 Fc q {v), with d v F^{v) = -F R q {v). One finds that 
at equilibrium and to this order in r, C and C coincide. 

The non equilibrium connected correlation is already 
non trivial in the absence of disorder I44L It also 



lim 



1 

xK 



= 2 + 2e? E T + 0(t 2 ) = 



(20) 



independently of v and e, i.e. of (small) wavevector and 
initial condition. There is thus a non trivial asymptotic 
FDR in the CO glass phase, which similarly to the case 
of pure critical points is in the interval [1/2, 0[. Here 
however it continuously de pen ds on temperature T. 

As in previous works |r3. lift ll2T | one can also examine 
the "diffusive" mode q = 0. It is possible to obtain a 
simple analytical form for any u — t/t' in the case e = 
(flat initial conditions <f>(x, t = 0) = 0), up to order 0(t 2 ) 



3q = ff ff («) = 2 + 2re 7E (l - atanh l-^j J (21) 



Although this quantity depends on e in general it reaches 
an e-independent limit for u — > oo. It is also interesting 
to compute X^ t , in real space. In previous works [6j, it 
was argued that 



lim lim Xfij 

t' — >oc t — *oo 



x=0 



lim lim X-> t 7° = Xoo (22) 

t' — ^oo t — >oo 
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The ^-independence found above (|20(l puts in the present 
case these heuristic arguments on firmer grounds 0. 

Following the discussion in Q , this Xoo leads to an ef- 
fective temperature T c g — T/X^, that can be measured 
by a thermometer coupled to the field 4>(x,t). Indeed, 
Fourier transforming Eq. (|12l I16fl . one checks that the 
local 7^7° and Cf t T° are precisely of the form given in 

It turns out that it is crucial to consider the connected 
correlation to obtain a non trivial FDR. We have also 
performed the calculation j^] for the correlation func- 
tion C\ v It exhibits a scaling form similar to i|ltj|) 
which decreases more slowly that C^ t , for large u. If we 
impose limu-toe uFc(v, u) — Fboc (v) one finds , for e = 
that C = 6 C + 1 , leading to A c = 1 - e< E r + O (t 2 ) . The 
FDR is found to be X^T = (2 + re 1E V") -1 and thus 
does not seem to approach (uniformly in r) a non trivial 
limit at large u. Thus, although both correlations give 
the same equilibrium result to order O(r), only the con- 
nected one, as defined here, yields a non trivial asymp- 
totic FDR in the non equilibrium regime. For e ^ 0, the 
large u behavior is dominated by the initial condition 
and Fc(v,u) ~ -Fboo(v) [42l lifij. The present consider- 
ations are also of interest for pure systems when initial 
conditions are non zero, e.g. in a quench from an ordered 
phase (or correlated initial conditions e ^ 0) 0, 142I lijj . 
Indeed, in that case one can distinguish connected and 
non-connected correlations (@J which can also lead to dif- 
ferent behaviours of the FDR. 

We have found a physically relevant glass phase in 
which one can show analytically the existence of a non 
trivial FDR. It is a robust quantity independent of the 
initial condition under study and, for some observables 
appear to be related to an effective temperature T e g. Its 
continuous dependence on T reflects the marginal char- 
acter of this glass phase. Our analytical predictions can 
be tested in numerical simulations [4p| . They are in good 
agreement near T g and it would be interesting to inves- 
tigate the aging behaviour at lower temperature, where 
hints of some new physics already appear in the statics. 
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